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ABSTRACT 

We calculate the orbital angular momentum of dark matter subhaloes in the Aquarius 
simulations of cold dark matter (CDM) galactic haloes. We calculate the orientation 
of their angular momentum relative to that of the spin vector of their host halo and 
find a variety of different configurations. All six Aquarius haloes contain statistically 
significant populations of subhalo orbits that are aligned with the main halo spin. 
All haloes possess a population of subhaloes that rotates in the same direction as 
the main halo and three of them possess, in addition, a population that rotates in 
the opposite direction. These configurations arise from the filamentary accretion of 
subhaloes. Quasi-planar distributions of coherently rotating satellites, such as those 
inferred in the Milky Way and other galaxies, arise naturally in simulations of a ACDM 
universe. 

Key words: galaxies: formation - dark matter 



(N 
> 

oo 

q 

od 
o 
o 



1 INTRODUCTION 

It has been known for several decades that the 11 'clas- 
sical' satellites of the Milky Way defin e a thin plane 
around the Galaxy l|Lvnden-Belll Il976l ). Some of the 
faint satellites recently discovered in t h e Sloan Digi- 
tal Sk y Survey (SPSS; lYork et all |2000|; IWillman et all 
2005a, b; Sakamoto & Hascgawa 2006; Zuckcr et al. 2006a b; 
Irwin et al.l 120071 : 1 Walsh et al.l 120071 ; iBelokurov et al.ll2008l ) 



also appear to have an ani sotropic distribut ion reminiscent 
of that of the classical 11 (|Metz et al.| [2009). The presence 
of such a 'disc-of-satellites' suggests a common plane of ro- 
tation in the Milky Way. Measurements of proper motions, 
which are now possible for some of the satellites, can be used 
to constrain the natur e of any systemic rotation (|Metz et al.l 
120081 ; iLux et alj|2010h . 

In tandem with these observational developments, ad- 
vances in computational cosmology now make it possible to 
simulate galactic haloes with sufficient resolution to probe 
the properties of satellites and investigate the origin of their 
flattened configuration. N-body simulations from cold dark 
matter (CDM) initial conditions show that a large number 
of accreted haloes survive to the present, making up a popu- 
lation of 'subhaloes' of the 'main halo,' some of which could 
host the satellites. 

The observations, however, suggest a complex forma- 
tion history. First, the number of satellites identified so far is 
much smaller than the number of dark subhaloes in the sim- 
ulatio ns, giving rise to the so-called 'missing satellite prob- 
lem' l|Moore et all 1 19991 ; iKtypin et al.| [l999l. Secondly, the 



thin ring around the sky delineated by the classical satellites 
contrasts with the dist ribution of subhaloes in the simula- 
tions which is triaxial (|Libeskind et al. I l2005l ; I Zentner et al.l 
2005). Thirdly, the inferred angular momentum vectors of 
the majority of the classical satellites in the Milky Way 
point towards a patch o n the sky of diameter n o greater 
than 30°, which has led iMetz et ail (|2008l . 120091 ) to argue 
that the observed satellites cannot have formed in cold dark 
matter subhaloes. In contrast to the Milky Way, NGC 5084 
appears to have a populatio n of satellites orbitin g in the 
opposite sense to the galaxy (Carignan et al. 1993)- 

The combination of 'missing satellites', an anisotropic 
distribution and coherent orbits is so metimes viewed as 
a challenge to t he CDM model (e.g. iMoore et all Il999l ; 
IMetz et al.l 120081 ) . However, a number of studies using 
semi-analytic modelling and hydrodynamic simulations have 
shown that a relatively small satellite population is a natu- 
ral outcome of ga l axy fo r mation in the CDM c osmology (e.g 
Kauffma rmet al.l [l993; Bu llock et al. | |2000t iBenson et al 



120021; ISomervilld 120021; iKoposov et al 



2009; 



Busha et alJ|2010t ICooper et alJ^OlOt iLi et alj|20ict 



l2Q09t iMunoz et al 



iMaccio et al.ll20101 ; IWadepuhl fc Springelll201ll ). The simu 
lations show that satellites form only in a small fraction of 
subhaloes which turn out to be those that h ad the most mas- 
sive progenitors at the time of accretion (fLibcskin d et al.l 
2005). Furthermore, disc-like subhalo configurati ons are seen 



to form in ACDM cosmological simulations ( Kang et al 



2005; Zentner et al. 2005; Libcski nd et~aiT 2007; 



Li fc Helmi 



2008; [Libeskind et al.l 20091 ). Such systems appear to be re- 
lated to the preferential accretion of haloes along the fil- 
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aments of the cosmic web. Haloes tend to fall along the 
central spines of filaments, so that the range of trajectories, 
and thus o rbits, that they acquir e when they enter a halo is 



rest ricted llLibeskind et al 



Shaw et al.l 



Libeskind et al.l 



2009). 



( 20061) . IWarnick fc Knebd |2006h and 
2009 ) confirmed the conclusion of 
Libeskind et al.l l|2005l ) that satellite accretion is a highly 
anisotropic process and found in their simulations a signifi- 
cant population of sub haloes that co-rota ted with the spin 
of their hosts. However, IShaw et al.1 (120061) simulated galaxy 
cluster haloes, not galactic haloes; Warnick Sz Knebel (2006) 
also focused on cluster haloes except for one example of a 
galaxy halo which, however, had only moderate resolution 
(a minimum subhalo mass of ^ min — 5.7 x 1O 7 M ). The 
largest sample of gala xy halo simulations so far is that of 
iLibeskind et al.l (|2009l ). They analysed 436 haloes but were 
only able to resolve subhaloes of mass m m i n ^ 2.76 x 1O 9 M0. 

In this study, we analyse the state-of-the-art, high res- 
olution simulations of six galactic haloes of mass ~ 1 x 
10 12 M o of the Aquarius project |Springel et al.1 l2008bh . 
These simulations resolve subhaloes of mass exceeding ~ 
1O 5 M0. We calculate the angular momentum of subhaloes, 
and use the results to interpret the Milky Way data. The pa- 
per is organised as follows. In Section [2] we briefly describe 
the Aquarius project and the analysis performed for this pa- 
per. The results follow in Section [3] and our conclusions in 
Section [4] 
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2 THE SIMULATIONS 

The Aquarius project is a set of dark matter simulations 
containing haloes similar in size and environment to those 
of the Milky Way; each one has been run from z = 127 to 
2 = 0. There are six different haloes (Aq-A - Aq-F), each of 
which has been resimulated at at least two levels of resolu- 
tion (LI, the highest, down to L5, t he lowest). They wer e 
performed using the GADGET- 3 code |Springel et al.ll2008bl ). 
In all cases, the resimulations at different resolutions show 
remarkable convergence in the positions and internal prop- 
erties of subhaloes. This project has already yielded several 
interesting results, incl uding a study of the near-universality 
of ha lo density profiles (|Springel et al.ll2008bl ; iNavarro et al.l 
2010), predictions for the 7 -ray signal from annih ilating dark 
matter in the galactic halo (jSpringel et a l. 2008a) and for di- 
rect d ark matter detection experiments ( Vogelsberger et al.l 
2009). A summary of key parameter values for each simula- 
tion is given in Table [T] 

It is important for this study to establish that the sam- 
ple of six Aquarius haloes can be considered at least ap- 
proximately representative of the population of Milky- Way 
mass haloes as a whole. The Aquarius haloes are all drawn 
from the same parent cosmological simulation, and it is pos- 
sible to address this issue directly for several properties. 
The spins, concentrations and formation histories of the 
Aquarius haloes are compar ed to the parent population in 
iBovlan-Kolchin et al.l (|2010T l. Broadly speaking, the proper- 
ties of the Aquarius haloes span the expected range for the 
population as a whole. We give the values of the halo spin, 
concentration and formation redshift, defined as the redshift 
when half the halo mass is assembled, in Table [T] Also in 
the table we list the shape axis ratios for the haloes, approx- 



Figure 1. The distribution of cos <9h-S f° r the six Aquarius haloes 
at resolution level L2, where £?h-S is the angle between the main 
halo spin and subhalo orbit vectors. The six are separated into two 
panels according to whether or not they exhibit an antiparallcl 
tail greater than 0.5. Top panel: results for Aq-A2 (orange), Aq- 
D2 (green), and Aq-E2 (light blue). Bottom panel: as above, but 
for Aq-B2 (blue), Aq-C2 (red), and Aq-F2 (purple). In both cases, 
the dashed line corresponds to an isotropic distribution. 



imating them as ellipsoids. The axis ratios are taken from 
Vera-Ciro et al. (in preparation) and calculated for ellipsoids 
which are determine d by applying the iterative method of 
lAllgood et all (2006) to the haloes with the substructure re- 
moved (actually to the 'main halo', defined below). The six 
haloes show a range of shapes and are typ ical for ACDM 
haloes (|Allgood et al.ll2006l ; lBett et al.ll2007l ). 

The halo membership of each particle is determined 
using the friends-of-friends (FOF) algorithm (|Davis et al.l 
Il985l ). The particles in each FOF group are, in turn, as- 
signed to self-boun d structures using the subfind code 
(|Springel et al.ll200ll ) . We call the largest of these self-bound 
substructures the main halo, and the remainder we call sub- 
haloes. A small proportion (< 1 per cent) of the particles 
within the FOF group are found to form a 'fuzz' that is 
not gravitationally bound to any other object; they are not 
considered any further. 

Our primary aim is to determine the orientations of 
dark matter subhalo orbits in the Aquarius simulations and 
compare the results with data for galactic satellites. We cal- 
culate the 'main halo spin', defined as the sum of the angular 
momenta of all main halo particles about their centre-of- 
mass. For each subhalo, we calculate the 'subhalo orbital 
spin', defined as the vector associated with the angular mo- 
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Name 


rn p [Mq] 


^200 [kpc] 


M200 [M©] 


n B 


A 


C NFW 


^form 


g 


P 


Aq-Al 


1.712 x 10 3 


245.76 


2.523 x 10 12 


197484 




16.11 


1.93 






Aq-A2 


1.370 x 10 4 


245.88 


2.524 x 10 12 


30177 


0.027 


16.19 


1.93 


0.866 


0.687 


Aq-A3 


4.911 x 10 4 


245.64 


2.524 x 10 12 


9489 




16.35 


1.93 


0.862 


0.688 


Aq-A4 


3.929 x 10 5 


245.70 


2.524 x 10 12 


1411 




16.21 


1.93 


0.844 


0.700 


Aq-A5 


3.143 x 10 6 


246.37 


2.541 x 10 12 


246 




16.04 


1.93 


0.830 


0.685 



Aq-B2 


6.447 X 10 3 


187.70 


1.045 x 10 12 


31050 


0.022 


9.72 


1.39 


0.820* 


0.839* 


Aq-C2 


1.399 X 10 4 


242.82 


2.248 x 10 12 


24628 


0.020 


15.21 


2.23 


0.711* 


0.770* 


Aq-D2 


1.397 x 10 4 


242.85 


2.519 x 10 12 


36006 


0.012 


9.37 


1.51 


0.846* 


0.901* 


Aq-E2 


9.593 x 10 3 


212.28 


1.548 x 10 12 


30372 


0.017 


8.26 


2.26 


0.898* 


0.674* 


Aq-F2 


6.776 x 10 3 


209.21 


1.517 x 10 12 


35041 


0.050 


9.82 


0.55 


0.700f 


0.866f 



Table 1. Selected parameters of the Aquarius simulations used in this paper. The simulation name encodes the halo label (Aq-A, B, 
and so on) and the numerical resolution level (1 to 5, hereafter LI, L2, L3, L4, L5). m p is the particle mass, r200 the radius of the sphere 
of density 200 times the crit ical density, M200 the halo m ass within r200i n s the number of subhaloes within the main halo, A the spin 
parameter as determined by Boylan-Kolchin et al. and q, p the halo shape axis ratios b/ a and c/b respectively (Vera-Ciro et al., 

in preparation). The axes are denned as a ^ b ^ c for ellipsoids determined using the method of lAllgood et all l l2006h . Values with * or 
f superscripts were calculated for haloes at resolution levels L4 or L3 respectively. As the smallest subhaloes determined by SUBFIND 
contain 20 particles, the minimum subhalo mass in each simulation is 20m p . 



2.0 




Figure 2. The distribution of cos 6>h-s f° r an the resolution levels 

of Aq-A. Aq-Al contains 197484 subhaloes (purple), Aq-A2 30177 fl f~l 

(orange), Aq-A3 9489 (red), Aq-A4 1411 (green), and Aq-A5 246 ' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 

(light blue). The error bars denote Poisson uncertainties. -1.0 -0.5 0.0 0.5 1.0 

cos H . S 



mentum of each subhalo about the centre of the main halo. 
We then calculate the cosine of the angle, Bu s, between the 
main halo spin vector and the subhalo orbit vector for every 
subhalo associated with that main halo. These subhaloes are 
tracked back to the initial conditions in order to investigate 
the origin of the patterns that we find. 



Figure 3. The distribution of cos 0h s f° r different populations. 
Top panel: comparison of the 1000 most massive subhaloes at 
z = (purple) with the 100 subhaloes that have the most mas- 
sive progenitors (light blue), and the entire population of Fig. [2] 
(orange). Bottom panel: comparison of the cos Sjj.s distribution 
for subhaloes with two populations of main halo particles: a sam- 
ple of 3 x 10 4 selected to have the same radial distribution as the 
main halo (blue) and the full population (red). 
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3 RESULTS 

We first describe our calculation of the angular momen- 
tum distributions of various populations and then investi- 
gate their origin. 

3.1 Angular momentum distributions of subhaloes 

We compute cos 8h s for each of the six L2 haloes at z — 
as described above, and plot the results in Fig.Q]as a proba- 
bility density; an isotropic distribution of angular momenta 
in this plot corresponds to a horizontal line at 0.5. 

All the haloes show a statistically significant bias for 
subhalo orbits to be aligned (para llel) to the rotati on of 
the parent main halo, as found by IShaw et al, I (I2006T ) and 
IWarnick fc Knebei §006). The average fraction of corotat- 
ing subhaloes in the Aquarius haloes is 57 per cent, with a 
narrow range between 54 per cent and 61 per cent. This 
is consistent with the a verage of 59 per cent quoted by 
IWarnick fc Knebei (|2006[K_ This result is a natural outcome 
of tidal torque theory (|Hovlelll95ll ; IWhitei H984T ) when the 
primordial dark matter protohaloes exert torques on one an- 
other, inducing net spins as they condense. 

We also find significant numbers of nearly anti-parallel 
orbits in three of our haloes. Specifically, haloes Aq-A2, Aq- 
D2, and Aq-E2 show a significant proportion of subhalo or- 
bits in the —1.0 <cos 8h s < —0.9 bin (9.5 per cent, 6.3 
per cent, and 7.3 per cent respectively where 5 per cent 
would be expected for a random distribution - the Pois- 
son errors on our L2 measurements are negligible), whilst 
Aq-B2, Aq-C2, and Aq-F2 do not. We have separated the 
haloes into two panels according to this property. We find 
an anti parallel excess in three o ut of six of our haloes, 
whereas IWarnick fc Knebei (2006) only have one such halo 
out of their sampl e of ni ne. Adopting the same binning as 
IWarnick fc Kne bc (2006) does not change our result. With 
such small halo samples it is unclear whether this particular 
result is consistent or inconsistent between the two studies. 

To test if our results are robust to changes in resolution, 
we repeat this calculation for the five different resolution lev- 
els of the Aq-A halo (Fig. [2}. We see that Aq-Al together 
with Aq-A3, Aq-A4, and Aq-A5 has an angular momen- 
tum distribution broadly of the same form as Aq-A2, with 
increasing noise as the resolution decreases because of the 
smaller number of subhaloes. Each resolution level is dom- 
inated by a different subhalo mass; the minimum subhalo 
mass in Aq-A5 is ~ 10 7 Mq, while in Aq-Al it is three or- 
ders of magnitude smaller. We find a similar degree of con- 
vergence with numerical resolution for haloes Aq-B through 
to Aq-F. 

In Fig. [3] we probe the orientation of the angular mo- 
mentum vector of different populations. In the top panel, 
we compare the distribution for the 1000 largest subhaloes 
at the final redshift (particle number > 1222, equivalent to 
subhalo mass of 1.7 x 10 7 Mq) with that the 100 subhaloes 
present at z = that had the most massive progenitors and 
that of the entire halo population. The most massive pro- 
genitor is defined as the SUBFIND halo in the merger tree 
that contained the largest number of particles over the en- 
tire history of the simulation. This mass is very close to 
the mass that the subhalo had at the time it fell into the 
main halo. It is these subhaloes that are most likely to host 



satellite galaxies, according to iLibeskind et~ai1 i|2009l ). Of 
the subhaloes that had the 100 largest progenitors, all bar 6 
are among the top 1000 most massive subhaloes at redshift 
zero. The distributions of cos #hs for all three populations 
of subhaloes are consistent within the errors. 

To establish whether the angular momentum orienta- 
tion of the subhalo population is special, in the lower panel 
of Fig. [3]we compare subhaloes in Aq- A2 with particles from 
the main halo. We create a special sample of halo parti- 
cles with the same radial distribution as the subhaloes. This 
is made by first defining a set of about 30 radial bins be- 
tween the halo centre and the virial radius. The halo sub- 
sample is produced by first noting how many subhaloes lie 
in a particular bin, and then randomly selecting the same 
number of halo particles from the that same bin. This is 
always possible as the number of halo particles in any bin 
exceeds the corresponding number of subhaloes. We com- 
pare this particle sample's distribution of cos #hs with that 
for the Aq-A2 subhaloes and for the entire set of main halo 
particles. The three distributions are statistically inconsis- 
tent with each other. The subhalo population has a larger 
fraction of aligned and antialigned members, with the radi- 
ally selected subsample being intermediate between the sub- 
haloes and the halo particles as a whole. Although even the 
latter has a non-uniform distribution of angular momenta 
cosines, it is significantly flatter than that of other two pop- 
ulations. This suggests that the accretion mechanism that 
supplies subhaloes (of all masses) is somewhat different from 
the mechanism by which halo particles are accreted, or that 
the evolution of subhaloes differs from that of halo particles. 

To investigate the orientation of the orbital spins in 
more detail, we plot the angular momentum vectors of each 
subhalo on an all-sky Mollweide projection, one for each halo 
at resolution L2. Each map displayed here was divided into 
~ 45000 pixels, with angular width ~ 1°, and smoothed 
with a Gaussian beam o f FWHM 10° using Healpix rou- 
tines (|G6rski et al.l |2005j) . We identify the pixel with the 
highest density after smoothing, and call this the 'densest 
point vector'. The pre-smoothing maps for all six L2 haloes 
are displayed in Fig [4] The main halo spin vector is marked 
in red, its antipole in blue, and the densest point vector in 
green. 

Aq-A2 exhibits the cleanest structure of all the haloes, 
with strong clustering around the pole and antipole, joined 
by two strands. Aq-B2 is, in contrast, characterised by ir- 
regular structures concentrated around regions distant from 
the main halo spin poles. All of the other haloes exhibit clus- 
tering around the main halo spin, with other, local, features 
apparent. The densest point vector position is always closer 
to the main halo spin than to its antipole. One may think of 
Figs.[T]toE]as an integration around lines of equal angle from 
the red and blue circles. As noted above, we are particularly 
interested in those subhaloes that are most likely to host 
satellites, and so we repeat this plot for the 100 subhaloes 
with largest progenitors in Fig. [5] 

As expected from Fig. [3] the 100 subhaloes with the 
largest progenitors trace the underlying structure of sub- 
haloes in the map traced in Fig. [4] A few of them lie in 
regions where there are few subhaloes of any mass, and so 
we might expect to find satellite galaxies spatially removed 
from the disc-of-satellites for at least some portions of their 
orbits. The majority, however, lie within underlying struc- 



© 2010 RAS, MNRAS 000, ??-Tlu1 



Linking Subhalo Accretion and Satellite Orbits 5 




Figure 4. Mollweide projections of the directions of the angular momentum vectors of all subhaloes in the L2 simulations. The red circle 
shows the direction of the main halo spin, blue the main halo spin antipole, and green the densest collection of vectors after smoothing. 
The maps have been rotated such that all three circles lie on the equator, with the main halo spin and its antipole lying 90° either side 
of the centre and the green circle in between. Thus, a subhalo of #hs = 0° w iH map to the red circle, and one of 0h s = 90° to either the 
plot boundary or a point on the north-south bisector. 



tures. The subhaloes with the top 11 progenitors cluster in 
the same way as the rest of the top 100. Thus, we conclude 
that observed satellite galaxies should also exhibit coherent 
motion. 



3.2 The origin of coherent rotation 

The importance of filamentary accretion can be appreci- 
ated by examining the positions of the subhaloes at different 
snapshots in the simulation. In Fig. [6] we plot the positions 
of all the subhaloes present at z — 0, relative to the centre 
of the main halo in two projections. On the left, the main 
halo angular momentum vector points along the positive X- 



© 2010 RAS, MNRAS 000, ??-fI01 



6 M. R. Lovell, V. R. Eke, C. S. Frenk, and A. Jenkins 




Figure 5. Mollweide projections of the directions of the angular momentum vectors of subhaloes with the largest progenitors in the L2 
simulations. Subhaloes with top 100 progenitors are denoted in blue; the subset with the 11 largest progenitors are plotted in purple. 



axis, so that the subhalo populations that have cos 8h s> 0.9 
(red) and cos 8n s< —0.9 (blue) appear as an edge-on thick 
disc. On the right, the angular momentum vector points out 
of the plane of the page. 

In Fig. [7] we investigate the origins of the different pop- 
ulations of subhaloes by plotting their positions in the initial 
conditions. No subhaloes have condensed at this early time, 
so we define the 'position' of each subhalo as the centre-of- 
mass of all the particles that will be members of that subhalo 
at redshift zero. Plotting the position of the most-bound 



particle rather than the centre-of-mass makes no difference 
to the appearance of the plot, and the plotting procedure 
followed is exactly the same as that used for Fig. [6] 

All of the haloes that have an excess of near-antiparallel 
subhaloes in Fig. [1] show a delineation in the positions of the 
different subhalo populations. There is also a degree of seg- 
regation in Aq-C2, however we find no clear delineation for 
haloes Aq-B2 and Aq-F2. We can observe how these segre- 
gated initial positions evolve into the orbital configurations 
at the final time by examining snapshots of intermediate 
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Figure 6. X — Y positions for subhaloes at z = in two projections. Subhaloes with cos 9n s> 0.9 are indicated in red and those with 
cos #h.s< — 0-9 m blue. All other subhaloes are shown in black. The black circle marks the virial radius, r200- Left panels (i.e. L): the 
X axis points in the direction of the main halo spin, so those subhaloes with orbit vectors parallel and antiparallel to the main halo spin 
appear as a band parallel to the Y axis. Right Panel (R): looking down the X axis, so the main halo is spinning anticlockwise. The red 
and blue points are plotted in a random order on top of the black. 



redshift. In the cases of haloes Aq-A and Aq-D, we find 
that, independent of resolution, the motion of interest of 
these delineated subhaloes occurs within the plane of the 
main halo spin (the panels A-R and D-R). This enables us 
to describe this motion simply with the schematic arrows 
shown in these panels in Fig. [Jj In both of these haloes the 
near-parallel and near-antiparallel populations collapse to 
form filaments in their segregated regions. Crucially these 
filaments are not straight, and the subhaloes follow curved 
paths into the main halo as shown by their same-colour ar- 
rows. 'Red' (near-parallel) subhaloes will therefore enter the 
virial radius with an anti-clockwise orbit around the plot 
centre, whereas the 'blue' will adopt a clockwise orbit. In 
Aq-D the 'red' and 'blue' filaments are separate entities, 
but in Aq-A they lie very close together and give the appear- 
ance of one filament fed at each end by two 'strands'. One of 
these strands then supplies the near-parallel subhaloes and 
the other the near-antiparallel. 



Aq-E near-antiparallel subhaloes are also accreted 
through a pair of curved filaments approximately in the 
plane of the final main halo spin, and so we illustrate the 
motion of these subhaloes with arrows in Fig. [7] panel E-R. 
By contrast, the accretion of the red subhaloes is more com- 
plex and involves motion at a significant angle to the plane 
of the main halo spin, and for this reason we do not draw 
the corresponding red arrows. In Aq-C some of the 'red' sub- 
haloes do accrete in a filament, but a large proportion end 
up in the large lump visible at the top of Fig. panels C- 
L and C-R. No coherent inflow pattern is apparent for the 
small population of 'blue' subhaloes. 

We can describe the accretion geometry further by de- 
termining where each subhalo enters the main halo. We find 
the redshift at which each subhalo attains its highest mass 
(taken to indicate the time when it falls into the virial ra- 
dius of the main halo) and thus determine its infall position 
relative to the main halo centre. The results are plotted in 
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Figure 7. X — Y comoving positions at z = 127 of the centre of mass of the particles that end up in each subhalo at z = 0. The 
coordinates are as in Fig. \E\ with the main halo spin at z = still determining the projections. The final virial radius is indicated in 
yellow. As stated in the text, the haloes Aq-A2 and Aq-D2 are found to accrete their 'red' and 'blue' subhaloes in the plane of the 
main halo spin, and so we have added coloured arrows to the A-R and D-R, panels to illustrate schematically the accretion paths for the 
different subhalo populations. The 'blue' subhaloes of Aq-E2 exhibit similar behaviour, and so we have also included arrows to indicate 
their motion in E-R. 



Fig. H3 which is oriented such that the main halo spin points 
towards the north pole of each projection. We can see that, 
independent of resolution, the subhalo populations that end 
up in parallel and antiparallel spin orbits in Aq-A, Aq-D, 
and Aq-E originate from preferential directions as expected 
from our visual examination. A majority of subhaloes in Aq- 
D and Aq-A accrete close to the equator, also as expected, 
whilst Aq-E acquires a significant proportion of its parallel 
orbit subhaloes from a patch of sky close to the main halo 
pole. Any demarcation for haloes Aq-B, Aq-C, and Aq-F is 
less clear, suggesting that filaments played a lesser role in 
their accretion history. 



4 CONCLUSIONS 

In this paper we have characterized the distribution of sub- 
halo orbits in the Aquarius simulations of CDM galactic 



haloes and attempted to explain the mechanisms that give 
rise to them. We find that the complex accretion patterns 
that build up a halo result in different configurations of 
subhalo orbits, none of which is close to isotropic. Some 
are structured in a symmetric way (Aq-A) relative to the 
spin poles, while others show no strong pattern (Aq-B). In 
all six haloes we find a large subhalo population that has 
coherent rotation aligned with th e spin of the ma i n halo , 
in agreement with the results of iLibeskind et al.l (2009). 
In three cases there is, in addition, a subhalo population 
that counter-rotates relative to the main halo. We trace 
this rather unexpected arrangement back to the filamen- 
tary nature of subhalo accretion. If gala xies tend to rotate 
in the same direction as their parent halo l|Bailin et al.|[2005l ; 
iBailin fc Steinmetzlkoorj ; iBett et aljkoid) , our results show 
that it is possible to generate populations of retrograde satel- 
lites. Such a population of retrograde satellites appears to 
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Figure 8. Mollwoide projections of the position at which each subhalo enters the main halo. Subhaloes that end up in parallel 
spin orbits are shown in red, those that end up in antiparallel spin orbits in blue, and those with intermediate orbits in black. 
The main halo spin points towards the north pole of each projection. Higher quality versions of these maps are available at: 
http:/ / astro, dur. ac. uk/^d50wse/page2. html 



be present in NGC 5084 |Carignan et al.lll997l ) whereas a 
population of pr ograde satellites a p pears to be pre s ent in 
the Milky Way (|Metz et alJl200St ); Hwang fc Parkl j201Ch 
find equal proportions of prograde and retrograde satellites 
across a sample of 215 systems. 

Halo Aq-A has a particularly concise formation history. 
This halo forms from a filament that collapses at early times 
and is fed by two strands at either end. A large fraction of 



the subhaloes that survive to the present pass though these 
strands, and are propelled into either a prograde or retro- 
grade orbit depending on the strand in which they origi- 
nated. Aq-D has a similar formation history and outcome, 
whereas Aq-E shows that it is possible to end up with a sim- 
ilar orbital arrangement by a different, more complex path. 
Coherent rotation is exhibited by the entire population of 
subhaloes, not just those with the most massive progenitors 
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which according, for example, to lLibeskind et ah! (|2005h are 
the most likely to host visible satellites. 

Our analysis has implications for the expected bulk 
kinematics of satellite galaxies which may be probed in fu- 
ture galaxy surveys. We expect a variety of orbital config- 
urations reflecting the variety of halo formation histories. 
Quasi-planar distributions of coherently rotating satellites 
should be commonplace, most rotating in the same direc- 
tion as the halo (and, by implication, the m ain galaxy) but 
some in the opposite direction as found bv lHwang fc Park! 
(l201(t) . 
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